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The first search for the rare kaon decay Kl — > has been performed by the E391a 

collaboration at the KEK 12-GeV proton synchrotron. An upper limit of 4.7 x 10 -5 at the 90 % 
confidence level was set for the branching ratio of the decay Kl — > 7r°7r°zyt> using about 10 % of the 
data collected during the first period of data taking. First limits for the decay mode Kl — > 7r°7r°P, 
where P is a pseudoscalar particle, were also set. 



PACS numbers: 13.20.Eb,11.30.Pb,12.15.Hh 

The decay Kl — > tPtPvv is a Flavor Changing Neutral 
Current process [l| involving a s — > dvv transition. In 
the standard model it is a predominately CP conserving 
mode with the branching ratio proportional to the square 
of the p parameter in the Wolfenstein parameterization 
of the CKM matrix. The predicted branching ratio in 
the standard model is (1.4 ±0.4) x 10~ 13 [2]. In this arti- 
cle, we report the first experimental limit on this mode. 
Additionally, we set limits on the supersymmetric mode 
Kl — > 7r°7r°P, where P is the pseudoscalar sgoldstino 

The spontaneous breaking of any global symmetry 
results in a massless Nambu-Goldstone mode with the 
same quantum numbers as the symmetry generator. In 
the case of supersymmetry, the symmetry generator is 
fermionic resulting in a Nambu-Goldstone fermion, the 
goldstino. The goldstino has a superpartner, the sgold- 
stino, which has scalar and pseudoscalar components. 
If the pseudoscalar sgoldstino is light enough (rap < 
mK L — 2ra 7r o) and the quark-sgoldstino coupling is parity 
conserving, then there should be the decay Kl — > 7r°7r°P. 
These conditions are met in a variety of models. An up- 
per bound on the branching ratio of this decay of « 10 -3 
can be derived from limits on the mass difference be- 



tween Kl and Ks [3[. Direct searches for the related 
charged modes K ± — » 7r ± 7r°P have set upper limits on 
those modes [4], [5]. The limit on the K + mode is 4x 10 -5 
for rap < 80MeV/c 2 . The K~ mode limit is 9 x 10~ 6 for 
rap < 200MeV/c 2 . The branching ratios of the charged 
modes are sensitive to the magnitude of the coupling be- 
tween quarks and sgoldstinos while Kl — > 7r°7r°P is sen- 
sitive to the real component of the coupling. Generally 
the K^ — > 7r ± 7r°P decay is suppressed by isospin conser- 
vation and therefore Kl — > 7r°7r°P has greater sensitivity 
to the coupling except when the imaginary component of 
the coupling dominates. 

The E391a experiment at the KEK 12-GeV proton 
synchrotron is a dedicated experiment for the decay 
Kl — > 7r uiy The first run of data taking took 

place from February to June 2004 (Run I). The anal- 
ysis in this Letter corresponds to approximately 10% of 
this data. The experiment uses a neutral beam extracted 
at 4° from the primary proton line. The beam was colli- 
mated into a circular beam with a 2 mrad half cone angle 
[8[. The mean K\ momentum is 3.5 GeV/c. 

The detector apparatus consists of a Csl crystal 
calorimeter and 4ir hermetic photon veto system. A di- 
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agram of the detector cross section is shown in Fig. [TJ 
The photon vetoes are arranged cylindrically around the 
beamline. The upstream end of the detector is 11 m 
downstream of the target, which we define as the origin 
of our Z -coordinate. 

The calorimeter consists of 576 blocks of undoped Csl 
crystals @. The majority of crystals are 7 x 7 x 30cm 3 
blocks which correspond to 16 X parallel to the beam 
line. There are 24 5 x 5 x 50cm 3 (=27X ) crystals sur- 
rounding the beam hole. The energy resolution was 
a E /E ~ 1%/y/E © 1% where E is in GeV, as measured 
with 25 Csl crystals and a positron beam. The aver- 
age position resolution of photon reconstruction is 5 mm. 
The face of the Csl array is located at a Z position of 
614.8 cm. In front of the Csl is a bank of scintillator 
plates (CV) for vetoing charged particles consisting of 32 
overlapping 6 mm thick plastic scintillator located at 550 
cm. 

There are multiple veto detectors making up the 47r 
photon veto system. The photon vetoes are the Front 
Barrel (FB), Main Barrel (MB), collar counters (CC02- 
07), and Back-Anti (BA) as shown in Fig. [U These are 
all lead-scintillator sampling calorimeters, except for the 
BA which also has layers of quartz. The Main Barrel is 
the primary photon veto surrounding the fiducial decay 
region. It is made of lead-scintillator sandwich with an 
inner diameter of 2.0 m and is 13.9 Xq thick. 

Most of the detectors are contained within a vacuum 
vessel. The fiducial decay region is kept at a pressure 
of 1 x 10 -5 Pa. This central region is separated from 
the rest of the detector by a 20 mg/cm 2 thick vacuum 
membrane. This membrane was not properly secured 
during Run I and hung into the beam line near the CV 
at ~ 50 cm upstream of the Csl face. The membrane 
acted as a target for neutron interactions. 

The photon vetoes were calibrated using either cosmic 
ray or muons from the target. The Csl array was cali- 
brated during a special run using a 5 mm Al target to 
produce 7r 0, s with a known position. This calibration is 
matched to cosmic ray muon data which is used to track 
changes in gain over time. 

The experimental signature of both Kl — > 7r°7r ^^ 
and Kl — > 7r°7r°P is four photons in the final state 
with missing mass. Their transverse momentum, Pt, 
has a higher maximum value compared to Kl — > 7r°7r° or 
Kl — > 7r 7r°7r° with missing photons. The energy and 
hit position of the four photons are measured by the 
calorimeter. Photons are reconstructed by summing the 
energy of contiguous blocks with energy deposited. The 
position of the photon is calculated by using the distri- 
bution of visible energy in the Csl blocks. 

The possible 7r° decay vertices can be reconstructed 
assuming the decay occurred on the beam axis and the 
7T° mass. There are multiple ways of pairing the photons 
to form 7r°'s, and for each pairing we calculate a x 2 of 
the difference in tt° Z vertices. We select the pairing with 



the minimum x 2 to determine the kaon Z decay vertex. 
The signal modes are distinguished by their relatively 
high Pt and reconstructed invariant masses of the 7r°-7r° 
system below that of the Kl- The Kl — > 7r°7r° mass 
peak is reconstructed with a = 8 MeV/c 2 . 

The kaon flux was calculated using Kl — > 7r°7r°. It was 
cross-checked by Kl — > 7r 7r°7r°. The acceptance of these 
modes, after applying selection cuts, was calculated using 
a GEANT3 based Monte Carlo simulation [10]. The sim- 
ulation included an overlay of accidental events selected 
from data. The acceptances and fluxes for these modes 
are shown in Table I. 

The signal region was defined by the Pt of the recon- 
structed kaon, the invariant mass of the 7r° — 7r° system, 
and the Z position of the reconstructed decay vertex. 
The acceptable Pt was defined to be between 100 and 200 
MeV/c. The lower bound was dictated by the presence 
of large amounts of 3tt° background at lower values of 
Pt . The invariant mass was required to be between 268 
and 450 MeV/c 2 . The lower bound is set by the minimal 
reconstructed mass with the intermediate reconstruction 
of two pions and the upper bound is set by the presence 
of Kl — > 7r°7r° mass peak at the kaon's true mass of 498 
MeV/c 2 . The acceptable decay vertices are between 300 
and 500 cm. The upstream limit of 300 cm is set by the 
presence of beam halo neutrons interacting with CC02 
at 275 cm. The downstream limit of 500 cm is set by 
core neutrons interacting with the vacuum membrane at 
~ 550 cm. 



Mode 


Acceptance 


Flux 


7T 7T 

_0 
7T 7T 7T 


(2.16 ±0.13) x 10" 4 
(1.39 ±0.07) x 10" 6 
(5.33 ±0.23) x 10" 5 


(1.54 ±0.04) x 10 y 
(1.57 ±0.04) x 10 9 
NA 



TABLE I: Acceptance and flux calculations of different 

signal modes. Acceptance is calculated from Monte 
Carlo. Flux is the number of kaon decays in the fiducial 
region. 

The acceptance of the mode Kl — > 7r°7r°P depends on 
the mass of the P. As the mass increases there is a decline 
in acceptance due to the reduction in the maximum Pt of 
the pions. This eventually leads to the phase space of 
the decay to lie completely in the region of high Kl — > 
7r 7r°7r° background. The single event sensitivity as a 
function of sgoldstino mass is shown in Fig. 2. 

There are two important sources of background to the 
signal: Kl — > 7r 7r°7r° with missing photons, and neutron 
related backgrounds. We predicted the backgrounds from 
data using a bifurcation technique U],[l2j]. ^he num ber 
of observed background events can be factored into the 
number of events after applying a group of setup cuts 
and the probability of one of these events surviving the 
remaining cuts. If we break up the remaining cuts into 
two uncorrelated groupings, A and B, this probability 
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FIG. 1: Cross section of the E391a detector. K^s enter from the left side. 
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FIG. 2: Single event sensitivity for Kl — > 7r°7r°P decay 
as a function of sgoldstino mass. 



can be factored into the probabilities of surviving each 
set, P(A) and P(B). We can rewrite this in terms of the 
number of events surviving the application of each set of 
cuts and the inverse of the other and the events which 
survive the application of the inverse of both cuts, 



Nb^ = N AB N AB /N AB . 



(1) 



7T 7T 7T . 



The primary source of background is Kl 
There are three mechanisms for this decay to produce 
background: two photons missing due to either geomet- 
ric or detector inefficiency, one photon missing and one 
photon fusion (two photons reconstructed as one) in the 
Csl, or two photon fusions in the Csl. Background events 
from all three mechanisms have similar distributions in 
the signal space. Our cut set A consists of photon veto 
cuts. Cut set B is made up of cuts on the quality of the 
photon cluster and the particle reconstruction. 

To check the bifurcation methodology, we applied it to 
regions surrounding the signal region. The Low Pt region 
is defined by same bounds in Z and invariant mass of the 
7r°-7r° system as the signal region and a Pt between 50 
and 100 MeV/c. The High Mass region is defined by 
the same bounds in Z and Pt as the signal region and a 
mass between 450 and 550 MeV/c 2 . The High and Low 
Z regions have the same bounds in invariant mass and 
Pt as the signal region and have reconstructed vertices 



between 500 to 550 cm and 250 to 300 cm, respectively. 
The predictions for these regions agree fairly well with 
data as shown in Table HI1 

For this technique to correctly predict the background, 
the cut sets A and B need to be uncorrelated. We se- 
lected cuts on this basis, but there is some correlation. 
An estimate of the uncertainty caused by ignoring this 
correlation can be made with values which are available 
without examining the signal region as 



C e =exN AB (l 



N 



)• 



(2) 



AB 



Here e is the difference in cut survival probability between 
cut A for an event passing cut B and an event passing the 
inverse of B. We determined e using events in the Low 
Pt region. We estimate this as a systematic uncertainty 
of 0.12 background events. 

One cause of cut correlation is contamination of other 
background sources in the signal region, primarily neu- 
tron related backgrounds. This is not measured by our 
estimate of e, because the neutron background is not 
present in the Low Mass, Low Pt region while it is 
present in the signal region. If there is a second back- 
ground source with different cut survival probabilities 
there is a correction to the prediction. It takes the form 
of 
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NabNab 
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(3) 



AB 



Here, N\ and N2 are the number of events of each back- 
ground type before cuts A and B are applied, and and 
Ab are the difference in the veto probabilities between 
the background types for the two cuts. It is important 
to note that the correction term does not directly corre- 
spond to the background contribution from the secondary 
sources. 

The largest source of neutrons is the interaction of 
beam core neutrons with the vacuum membrane in front 
of the CV. A second source is the interaction of halo 
neutrons with CC02. Both of the sources produce re- 
constructed events which are localized to their point of 
interaction with high Pt. Both of these regions are out- 
side the fiducial decay region. To estimate the impact 
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of these events on the the background prediction, it was 
necessary to determine the number of core neutron back- 
ground events in the signal region under the setup cuts. 

The number of neutron events in the region with Z 
greater than 500 cm and Pt greater than 0.1 GeV/c is 
too small to fit when all cuts are applied. We therefore fit 
the data using a loose set of cuts. We remove the BA cut, 
cuts on the distribution of photon energy and timing, and 
a veto on additional energy in the Csl. Additionally, we 
apply the inverses of cuts A and B to ensure we are not 
observing events in the signal region while being able to 
look at events in the Z fiducial region. With these sets 
of cuts the core neutron peak in the high Pt -high Z 
region is clearly visible and can be fit with a Gaussian 
as shown in Fig. 3. In the high Pt -fiducial Z region, 
there is a predominately Kl — > 7r 7r°7r° background. The 
Gaussian component of the core neutron peak produces a 
negligible contribution to the background in the fiducial 
region. The density of this distribution is calculated by 
subtracting off the Kl — > 7r 7r°7r° contribution found by 
Monte Carlo simulation and fitting the remainder by a 
Gaussian plus a straight line. Integrating this function 
over the fiducial region gives 248 ± 4 sta t. ± 124 syst . core 
neutron background events with the loose cuts applied. 
The application of the rest of the setup cuts reduces this 
by a factor of « 115. A predicted total core neutron back- 
ground in the signal region is 2.16 ± 0.03 s t a t. ± l-05 sys t. 
under the setup cuts, before the application of cuts A 
and B. This corresponds to N2 in Equation [3l 

The cut survival probabilities do differ significantly 
with A A = (26.4 ± 1.3)% and A B = (8.6 ± 0.5)% which 
are derived from Monte Carlo studies. Ni is estimated 
by assuming the background before the setup cuts are 
applied is dominated by Kl — > 7r 7r°7r° events and sub- 
tracting off the core neutron contribution. We estimate 
Ni = 102 ± 10. This results in an estimation of the un- 
certainty to the background prediction of 0.06 events. 

TABLE II: Prediction of background events in different 
regions. Statistical errors are shown. 



Region 


Nab 


N A b 


N AB 


Prediction 


Data 


Low Pt 


393 


72 


115 


21.1 ±3.3 


13 


High Mass 


46 


9 


4 


0.78 ± 0.48 


1 


Low Z 


5 














High Z 








6 








Signal 


84 


18 


2 


0.43 ±0.32 


1 



The bifurcation background predictions gives us a total 
background prediction of 0.43 ± 0.32 sta t. ± 0.13 sys .. In 
the signal region we observe a single event shown in Fig. 
H The event has M^o.^o = 315 ± 8 MeV/c 2 , P T = 
111 ± 5 MeV/c, and Z = 357 ± 11 cm. This is consistent 
with the background prediction. 
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FIG. 3: Z distribution of reconstructed events with 
Pt > 100 MeV/c. Points are data, shaded region 
represents a scaled Kl — > 7r 7r°7r° Monte Carlo sample, 
and the line is the fitted curve. The top plot shows the 
distribution with the loose set of cuts from which we 
derived our neutron background prediction, and the 
bottom plot shows the same distribution with all cuts 
applied. 



The uncertainty in the branching ratio of Kl — > 7r°7r° 
contributes a systematic uncertainty of 0.5% to our sin- 
gle event sensitivity. Uncertainty in the calibration of 
the photon vetoes contributes a systematic uncertainty of 
3.7%. Discrepancies between Monte Carlo and data ac- 
ceptance loss in the cuts, which were simulated by Monte 
Carlo, give a systematic error of 3.4% in acceptance. The 
total systematic acceptance error is 5.0%. Combining the 
branching ratio and acceptance errors produces a 7.1% 
systematic error in the single event sensitivity. 

Our single event sensitivity for Kl — > 7T 7r°ufy is 
(1.20±0.06 st at.±0.09 sys .) x 10~ 5 . With 1 observed event, 
which is consistent with the background prediction, we 
set a limit at the 90% confidence level for the branch- 
ing ratio at 4.7 x 10 -5 using Poisson statistics. This is 
the first limit set on this decay mode. The single event 
sensitivities for Kl — > 7t°7t°P with different sgoldstino 
masses are shown in Fig. O For mp < 100MeV/c 2 , we 
can set a limit on K L -> 7r°7r°P of 1.2 x 10~ 6 at the 90% 
confidence level. 

This analysis uses 10% of the data from the first E391a 
run. It appears that with this data set the analysis is 
background limited. We expect the second and third 
runs of E391a to have significantly reduced neutron back- 
grounds due to correcting the vacuum membrane issue 
and an insertion of a Be absorber into the beam. The 
second and third runs have more than 10 times the data 
as this analysis. However, further study is needed to de- 
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FIG. 4: Pt plotted versus invariant mass with all 
selection cuts applied. The rectangular regions 
correspond to the regions in Table HI The 2tt° signal 
region was used to calculate the kaon flux in Table I. 

termine the proper cut selection and what acceptance can 
be achieved in this improved neutron environment. 
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